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Abstract 

A mechanism for explosive magnetic reconnection is investigated by analyzing the nonlinear 
evolution of a collisionless tearing mode in a two-fluid model that includes the effects of electron 
inertia and temperature. These effects cooperatively enable a fast reconnection by forming an X- 
shaped current-vortex layer centered at the reconnection point. A high-resolution simulation of this 
model for an unprecedentedly small electron skin depth de and ion-sound gyroradius ps, satisfying 
de = Ps, shows an explosive tendency for nonlinear growth of the tearing mode, where it is newly 
found that the explosive widening of the X-shaped layer occurs locally around the reconnection 
point with the length of the X shape being shorter than the domain length and the wavelength of 
the linear tearing mode. The reason for the onset of this locally enhanced reconnection is explained 
theoretically by developing a novel nonlinear and nonequilibrium inner solution that models the 
local X-shaped layer, and then matching it to an outer solution that is approximated by a linear 
tearing eigenmode with a shorter wavelength than the domain length. This theoretical model 
proves that the local reconnection can release the magnetic energy more efficiently than the global 
one and the estimated scaling of the explosive growth rate agrees well with the simulation results. 
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I. INTRODUCTION 


Clarification of the mechanisms of fast magnetic reconnections in space and laboratory 
plasmas is a fundamental issue that has been tackled by plasma physicists for more than 
half a century HIE]. In particular, the explosive release of magnetic energy, observed 
in solar flares, magnetospheric substorms and tokamak sawtooth collapses, indicates that 
the magnetic reconnections might further accelerate through nonlinear and nonequilibrium 
processes. It is widely accepted that the resistive magnetohydrodynamic (MHD) model 
cannot reproduce such fast reconnection unless a locally enhanced resistivity is artihcially 
introduced [3]. Since the plasmas in these explosive events are relatively collisionless, the 
resistivity is thought to be physically less important than non-collisional microscopic effects 
such as electron inertia. Hall current, ion gyroradius effects, and so on, which are all neglected 
in the classical MHD model. Consequently, there has been considerable effort in recent years 
in studying a variety of two-fluid W\, gyro-fluid mi and gyrokinetic mm models to 
understand reconnection in collisionless plasmas. Until now, many simulation results 
flj have shown that collisionless reconnection tends to accelerate into a nonlinear phase. 
However, the theoretical understanding of this process is very limited. 

The primary computational obstacle is that the nonlinear acceleration phase is observed 
only when the magnetic island width (or the amplitude of the tearing mode) exceeds the 
microscopic scales while sufficiently smaller than the scale of the equilibrium magnetic shear. 
This suggests the importance of making the microscopic scales as small as possible, yet 
affording very high spatial resolution to prolong the acceleration phase. Consequently, it 
has not been clear in previous simulations how reconnection is accelerated and whether or 
not it is explosive. 

The nonlinear theory for explosive magnetic reconnection has remained elusive because 
of the difficulty of solving strongly nonlinear and nonequilibrium fluid motion with multiple 
scales. The method of asymptotic matching has been only applicable to linear stability of 
collisionless tearing modes, where the island width is assumed to be much smaller than any 
microscopic scale [HI 1201122] . In the dissipationless limit, recent studies take advantage of 
the Hamiltonian structure [26II28] of the collisionless two-fluid models [I11II81129H3I]. These 
studies show that the two-fluid effects distort the conservation laws (frozen-in variables) and 
hence permit magnetic reconnections with ideal fluid motion. 
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In the presence of only electron inertia in the two-fluid model, an ideal fluid motion 
develops an elongated current layer with Y-shaped ends, where the layer width is comparable 
to the electron skin depth (de) [13] • On the other hand, even faster reconnections due to 
the formations of X-shaped current-vortex layers are observed numerically naE] when the 
effect of electron temperature is taken into account and the ion-sound gyro-radius {ps) is 
comparable to or larger than de- This distinction between the Y and X shapes seems to 
be crucial in determining the reconnection speed, in analogy with that between the Sweet- 
Parker pQ and Petschek [2] models for resistive reconnections. Although several pioneering 
works [131 ES] have attempted theoretical explanations of the explosive growth of these 
nonlinear tearing modes, we note that their predictions are not in quantitative agreement 
with the high-resolution simulation results given in Ref. pa and the present work. 

The goal of this paper is to clarify an explosive mechanism for collisionless reconnection 
caused by the interplay of the effects of electron inertia and temperature. To this end, we 
analyze the simple Hamiltonian two-fluid model given explicitly in Sec. both numerically 
and analytically. In previous work [23] , we considered only the effect of electron inertia and 
estimated an explosive growth rate by using a new variational method. This method not 
only gives better agreement with simulation results than earlier work ra, but also gives a 
better physical interpretation because energy conservation is properly taken into account. 
Here we generalize our previous study for the Y-shaped layer and consider an X-shaped 
layer. 

To be more specihc, we will restrict our consideration to the case of de = Ps, for simplicity, 
and shorten the scale de = ps as much as possible in the simulations. For the hrst time we 
perform simulations with ps = dg < O.OIL and hnd that the X-shaped current-vortex layer 
widens rather locally around the reconnection point regardless of the size of computational 
domain. We show theoretically that this local X-shaped structure is nonlinearly generated 
because it is optimal for releasing the magnetic energy more efficiently than global structures. 

Our variational method [23] is inspired by the ideal MHD Lagrangian theory, in which the 
magnetic energy is considered to be part of the potential energy of the dynamical system 
(in analogy with the elastic energy of rubber bands). If a fluid displacement continually 
decreases the magnetic energy, it is likely to grow by gaining the corresponding amount of 
the kinetic energy, and the most unstable displacement would decrease the magnetic energy 
most effectively. This argument assumes that the two-fluid effects are essential for changing 
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the topology of magnetic field lines within the thin boundary layer, but their impact on the 
global energy balance is negligible in the limit of de, Ps —t 0. By choosing a fluid displacement 
as a test function that mimics the local reconnection, we can estimate a growth rate for the 
displacement from the kinetic energy. 

In Sec. 1^ we hrst introduce a reduced two-fluid model [5] that includes the effects of 
electron inertia and electron temperature, and focus on a collisionless tearing mode that is 
linearly unstable for a magnetic shear By{x) oc sin(27ra;/L3,) in a doubly periodic x-y plane, 
where the wavenumber ky = 2tt / Ly is related to the aspect ratio LyjLx of the domain. In 


Sec. Ill, we present our numerical results on the nonlinear evolution of this tearing mode. 
Explosive growth is observed when > 0.65 (where de = Ps and the tearing index A' is 
a function of Ly/L^). We will End that the explosive growth rate is almost independent of 
LyjLx because, in the explosive phase, the X-shaped current-vortex layer expands locally 
around the reconnection point and its characteristic length in the p-direction is shorter than 


Ly. In Sec. IV, we present our theoretical model that explains the generation of such a 
local X-shaped layer. By replacing Ly with a shorter length Ay(< Ly)^ we introduce an 
effective tearing index A' as a free parameter of the external solution and connect it to a 
novel inner solution that represents nonlinear evolution of a X-shaped layer. We will show 
that a local reconnection < Ly can transform the magnetic energy into the kinetic energy 
more efficiently than the global one Ly = Ly. Using this variational principle we estimate 
a growth rate for this local reconnection model, which is indeed explosive and agrees with 
the simulation results. We finally summarize in Sec. |V] 


II. MODEL EQUATIONS 

We begin with the reduced two-fluid model given in Refs. m with the assumption of 
cold ions, 

^ + [0,vV] + [vV,^] = o, (1) 

^ + [0 - 0] - dl + [0, = 0, (2) 

which governs the two-dimensional velocity field u = x \/(j){x,y,t) and magnetic field 
B = y/pQmino'V'ip{x, p, t) x + BqEz, where the guide field Bq and mass density rriino are 
assumed to be constant, po is the magnetic permeability, and [f,g] = (V/ x Vp) ■ is 
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the Poisson bracket. Here, the effects of electron inertia and electron temperature introduce 
two microscopic scales: the electron skin depth = cjupe and the ion-sound gyroradius 
Ps = \/%/ mi/ujci, respectively (where c is the speed of light, ujpe is the electron plasma 
frequency, Ud is the ion cyclotron frequency, Te is the electron temperature, rrii is the ion 
mass). 

The equations ([^ and (|^ conserve the total energy (or the Hamiltonian) that is given 
by the following energy integral, 

E = ^jd^x (IV0|2 + p^ivvr + IVV'P + dl\V^^l^\^) 

=: Ey + Et + Eb + Ec, (3) 


where Ey = f |V0p(i^x/2 is the ion perpendicular flow energy, Et = f p^l 
is the electron thermal energy, Eb = J \ V'ip\‘^d‘^x/2 is the magnetic energy and Ec = 
J d'l\W‘^'ip\^d?x/2 is the electron parallel flow (or current) energy. Kuvshinov et al. [29] 
and Cafaro et al. [Tl| show that ([^ and (|^ can be further rewritten as 


d'ip+ 

dt 

d'ljj- 


+ [0+,V’+] 
+ [0-,^-] 


0 , 

0 , 


(4) 

(5) 


in terms of 


ii± = 'tp - dlv‘^tp ± psdeV'^cj), (6) 

<P± = <p-plV^<P±PsdeV^ij. (7) 


It follows that and ip- are frozen-in variables, whereas ip is not unless de = 0. Magnetic 
reconnection is therefore possible when de 7^ 0 without any dissipation mechanism. 

As is common with earlier works miiisiiao], we consider a static equilibrium state, 

(p^^\x^y) = and ip^'^\x^y) = ipQCOs{ax)^ (8) 

on a doubly periodic domain D = [— La;/2, La,/2] x [—Ly/2,Ly/2] (where a = 2 tt/Lx), which 
is unstable with respect to double tearing modes whose reconnection layers are located at 
X = 0 and x = -ELx/2. For initial data we assume a sufficiently small perturbation of a 
single harmonic, (p oc sin kyy and ip — ip^^'^ oc coskyy, where ky = 21:1 Ly. Then, the following 
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parities, 


<Kx, y, t) = -iH-x, y, t) = -it,(x, -y, t), 
i>(x,y,t) = i>(-x,y,t) = tli(x,-y,t), 


(9) 

( 10 ) 


are exactly preserved by Eqs. 0 and ([^ for all t [32]. Therefore, the origin {x^y) = (0,0) 


(and the fonr corner points of D as well) is always a reconnection point. The solntions to 
this problem are fnlly characterized by three parameters; de/L^, Ps/Lx and the aspect ratio 

Ly/Lx{— cx/ky). 

The linear stability of this collisionless tearing mode has been analyzed in detail by many 
anthors For a given wavennmber ky in the y-direction, the tearing index at the 

reconnection layer a; = 0 is calcnlated as 



( 11 ) 


and the tearing mode is nnstable when A' > 0, namely, 0 < ky/a = L^jLy < 1. 

For Ps > de, the analytic dispersion relation [T8| predicts that the maximnm growth rate 
occnrs when 




( 12 ) 


Since de/L^ -C 1 and ps/L^ -C 1 are nsnally of interest, this is often very large. 

If it belongs to the range L^A' > 100 (or ky/a = L^/Ly < 0.377) in which A' is well 
approximated by 


L^A' ^ WLl/Ll = IQa^kl, 


(13) 


we can estimate the maximnm growth rate ymax at the wave nnmber ky^^s,x as follows. 


O^max 



(14) 


( 15 ) 


where 

Given this backgronnd material we now tnrn to our numerical simulations. 
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FIG. 1: Contours of ^l) and (j) when e = 5ps, where ps = de = 0.02Lx, LyjLx = vr. 

III. NUMERICAL RESULTS 

Equations Q and ([^ are solved numerically for various parameters using the spectral 
method in both the x and y directions and 4th-order Runge-Kutta for time evolution. The 
nonlinear acceleration phase is observed when the magnetic island width becomes larger than 
the reconnection layer width that is of order for the case Ps > de- To observe this 

phase for a longer period, we have performed all simulations with dg = Ps and narrowed the 
layer width (~ dg = Ps) as much as possible. The most demanding case dg = ps = O.OOSLa, 
requires 8192 x 8192 grid points in wavenumber space. 

The nonlinear evolution of Q and (|^ was studied in earlier works [m UHl [301132] , and 
we reproduce the main features, as shown in Figs. [^and[^ Since "0+ and '0- are frozen-in 
variables, their contours preserve topology as seen in Figs, [^a) and (b) (where the arrows 
depict the fluid motions generated by 0±). Then, spiky peeks of -0+ and '0_ are formed 
and their ridge lines look like the shapes of “\” and “/”, respectively, around the origin. In 
light of the dehnition (|^, the current and vorticity distributions can be directly calculated 
by using ^l)± and shown in Figs. |^a) and (b), which exhibit a “X”-shaped current-vortex 
layer dnii] whose width is of order dg = ps- We also note from Figs. |^c) and (d) that a 
relation dg\'V‘^'ijj\ ~ Ps\'V^(j)\ holds inside the layer (except at the reconnection points). 

As indicated in Fig. [^c), we denote the maximum displacement of the field lines at 
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(a) -V'V’ 



- 0.5 0 0.5 - 0.5 0 0.5 

xjLx xjL^ 


FIG. 2: Intensity distributions (red: positive, blue: negative) of the current vorticity 

and —± ps'V'^cj) when e = 5ps, where ps = de = 0.02La;, LyjLx = tt. 



FIG. 3: Parameters ky and Ps{= de) that result in explosive growth (x and *) and saturation (□). 
The maximum linear growth rate occurs at points * for each pg. 

X = ±La,/4 by e and numerically measure it from the displacement of the contour ijj = 0. 
We have run simulations with various combinations of LxjLy and Ps/Lx{= de/Lx), and 
investigated whether e grows explosively or not. Our results are summarized in Fig. 
where it should be recalled that the tearing mode is linearly unstable at all points in this 
hgure since 0 < LxjLy < 1. The linear growth rate achieves its numerical maximum at 
points indicated by the asterisk (*) for each hxed which agrees with the theoretical 
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FIG. 4: Logarithmic plots of the displacement e versus its time-derivative e, where de = Ps = 0.01, 
and the current distributions (a) at e = 0.002La; for ky/a = 0.477, (b) at e = 0.047^2, for ky/a = 
0.239 and (c) at e = O.OdL^, for ky/a = 0.119. 


prediction (12) [and also (14) for ky/a < 0.377]. The sqnare symbol (□) indicates points 
where exponential growth e oc (with the linear growth rate 7 ) stalls before e reaches ps- 
On the other hand, at the crosses (x) and asterisks (*) the exponential growth is accelerated 
when e gets larger than ps- These two regimes seem to be divided by a curve ~ 0.65. 

The above mentioned tendencies are demonstrated in Fig. which is a logarithmic plot 
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FIG. 5: Contours of (a) cf) and (b) V' at e = O.OlSLa; = 2>ps and (c) current distributions at 
e/ps = 1,2, 3,5, where de = Ps = O-OOSL^,, ky/a = ky^^^Ja = 0.171. 

of e versus e for the case Ps/Lx = 0.01. For ky/a = 0.477 and 0.557, which belong to the 
satnration regime < 0.65, the cnrrent-vortex layer spirals aronnd the 0-point as shown 
in Fig.U a) and the growth of e decelerates. Althongh this occnrs in an early nonlinear phase 
e < Ps in onr resnlts, we note that the satnration mechanism is similar to the one fonnd by 
Grasso et al. [32|, namely, the phase mixing of the Lagrangian (or frozen-in) invariants 'ip± 
leads to a new “macroscopic” stationary state. For the cases of ky/a = 0.119, 0.239, and 
0.398, which belong to p^A' > 0.65, we observe a transition from the exponential growth 
e cx e to an explosive growth e oc e" (n > 1 ) aronnd e ~ ps, and the latter continnes nntil 
the reconnection completes at e = L^/^. We farther note that, for the small ky/a = 0.119, 
a local X-shaped layer is generated spontaneonsly aronnd the reconnection point and it 
expands faster than the global one that originates from the linear eigenmode [see Fig. ic)]. 
By comparing the case ky = 0.119 with ky = 0.239 at the same amplitnde e = O.OdL^,, 
we hnd that this local X-shaped strnctnre aronnd the origin in Fig. |^c) is identical to 
that in Fig. [^b). Therefore, the explosive reconnection seems to be attribnted to the fast 
expansion of the local X-shape with a certain optimal size that is independent of the domain 
length Ly = ‘li^/ky. In fact, the nonlinear growth rates for ky = 0.119, 0.239, and 0.398 are 
eventnally comparable for e > 0.02 becanse the released magnetic energies are almost the 
same regardless of ky. 

We remark that the length of the local X-shape is not simply related to the wavelength 
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27r//cy,max of the most linearly unstable mode. For de = Ps = O.OOSLa; and ky = = 

0.171a, Figs, [^a) and (b) give a closer look at the contours of (j) and 'ifj, respectively, at 
e = 3ps, and Fig. |^c) shows the shapes of the current layer observed at e/ps = 1,2,3, and 5. 
As can be seen from Fig. |^c), a local X-shape appears and expands quickly in the nonlinear 
phase e > ps even though this reconnection is triggered by the most linearly unstable tearing 
mode ky = ky^^nax- Under the same conditions, the evolution of the energies Ev,t,b,c dehned 
by (|^ is shown in Fig. (where the total energy conservation E = Eq = const, is satished 
numerically with sufficient accuracy). In the linear phase e Ps = O.OOSLa,, the magnetic 
energy E^ is transformed into Ey, E^ and Ec at different but comparable rates. However, in 
the nonlinear phase e > ps, we note that the magnetic energy is exclusively transformed into 
the ion flow energy Ey and the energy balance 6Ey + 6Eb ~ 0 is satished approximately. 
Since 6Et and 6Ec are negligible, we infer that the nonlinear dynamics is dominantly 
governed by the ideal MHD equation (de, Ps —t 0). This fact motivates us to regard K := Ey 
and W := Eb as the kinetic and potential energies, respectively, according to the MHD 
Lagrangian theory. 


IV. THEORETICAL MODEL 


In this section, we develop a theoretical model to explain the explosive growth of e in the 
nonlinear phase dg = Ps e Lx/^- The current-vortex layers are obviously the boundary 
layers caused by the two-huid effects and their width should be of order de = Ps- The ideal 
MHD equations, ([^ and ([^ with dg = Ps = 0, would be satished approximately outside 
the boundary layers. Moreover, we also assume that 0 and 0 are continuous across the 
boundary layers, because Ec = f dg|V^0pd^a;/2 and E^ = f p^| V^0pd^a;/2 are negligible 
in the energy conservation (Fig. in the limit of de,Ps —t 0. Note that V'tp and V0 may 
be discontinuous across the layer because it is a current-vortex layer. 

Based on these assumptions, we consider a family of virtual displacements that generates 
a local X-shaped current-vortex layer, and then seek the displacement that decreases the 
magnetic (or potential) energy most ehectively. Our reconnection model is illustrated in 
Fig. 0 where we show only the hrst quadrant around the origin owing to the parity ( |Io| ). In 
Fig-ig the magnetic held lines are assumed to be piecewise-linear and the red line denotes 
the boundary layer (i.e., the upper right part of the X-shape). 
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FIG. 6: Changes of energies Ey, Eb, Eq and Et versus e, where = Ps = O.OOSLj;, ky/a 
^y,max/® — 0.171. 



FIG. 7: Local reconnection model. 
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We will mainly focus on the three regions: (i) boundary layer (ii) inner solution with a 
X-shaped layer (hi) external solution. These regions are characterized by four parameters 
Ay, cr, /a;, ly as follows. First, the position of the boundary layer is specified by ly). Second, 
the displacement of the held line that is about to reconnect at the origin is denoted by a, 
which can be also regarded as the half width of the local island. Finally, to allow for local 
reconnection, we introduce the “wavelength” Ay of the displacement at x = Ta;/4, which 
may be smaller than the wavelength Ly{= 27r/ ky) of the linear tearing mode. We will assume 
the following orderings among these parameters: 

de = Ps ^ < a L^/A^ lx < cr ly < d^y/A < Ly/A. (16) 


A. Matching conditions across the boundary layer 


First, we focus on a neighborhood of (i) the boundary layer and introduce a local coor¬ 
dinate system {X, Y) in the frame moving with the boundary layer, so that the X and Y 
directions are respectively normal and tangent to the layer (see Fig. [^. Let the inner region 
of the boundary layer be —5 < X < 5, where 5 ~ de = Ps- In this coordinate system, the 
velocity v and the Alfven velocity b = B/^/pqPq are assumed to be uniform outside the 
layer. Using the continuities of "0 and (j) across the layer, we assign linear functions. 


^ - bf'^x + 

= (j)c + v[‘^'^X - VnY, 

on the down-stream side {X < —6) and 

V = - hPx + h,y, 

= 4>, + vPX - VnY, 


(17) 


(18) 


on the up-stream side {X > 5), where all coefficients depend only on time. The disconti¬ 
nuities of the tangential components, iJf^ ^ and ^ indicate the presence of a 
current-vortex layer within [—5,(5]. 

Since we have assumed that dti^ + [0, '0] = 0 holds outside the layer, [0, 'ip] must be also 
continuous, namely. 


- Vnbf’ = “ Vnbf\ 


(19) 
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FIG. 8: Neighborhood of the current-vortex layer. 


is one of the matching conditions between (17) and (18). 


Moreover, the conservation laws of 'ijj± require that '0_ has a spiky peek within the 
boundary layer whereas does not [see Figs, [^a) and (b)]. Since -0+ ~ holds, we hnd 
a relation, 

— ps'^'^4> on [—5,5], (20) 

which specihes the ratio between a current peak (—and a vorticity peak (V^0) inside 
the layer. We have already noticed this relation in Fig. By integrating these current and 
vorticity distributions over [—5, 5], we obtain another matching condition. 


-Ps{vl 


^) = d,{hf-bf). 


( 21 ) 


When the boundary layer is moving at a speed W in the X direction, the condition (19) 
is transformed to 

vf^hn - {Vn - Vn)hf^ = V^'^^bn “ (Un “ 14)5;''“^ 




( 22 ) 


in the rest frame, and the condition (21) is unchanged. Since these conditons also yield 




(23) 


we need to impose at least two matching conditions among (21), (22) and (23). 


B. Modeling of the X-shaped boundary layer 

Next, we consider (ii) the inner solution that contains the X-shaped boundary layer. The 
detailed sketch of this region is given in Fig. where the displacement map (xq, yo) {x, y) 
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FIG. 9: First quadrant of the X-shaped boundary layer. 


on the up-stream side (i.e., the right side of the boundary layer) is simply modeled by 


a 


X = xq +-iyo-ly) and y = y^. 

t"j/ 


(24) 


This displacement map fully determines '0 and (j) on the up-stream side as follows. Since 


a Lx/4: is assumed in (16), we expand the equilibrium flux function, 


'ip^^\x,y) = 'ipocos{ax) = V’o ( 1 — ^x^ ) +0{x'^), 


a 


(26) 

and neglect 0{x^) in this region [0 ,ct]. Except on the boundary layer, the magnetic flux ijj 
is frozen into the displacement and hence becomes 


'4)^^\x,y,t) ='4)^^\xo,yo) 


a 


x-yiy-ly) 

til 


(26) 


on the up-stream side. By regarding the parameters <jit),lx(t),ly(t) as functions of time, 
the time derivative of the displacement map gives the stream function, 

,2 


= -s,|- +179, 


(27) 


where Sy = cr/ly^ and the parity 0(a;, 0,f) = 0 has been used as the boundary condition. 

Now, let us consider a magnetic held line that is labeled by its initial position x = xq 
(where < Xq < a). When this held line is displaced by the map, it intersects with the 
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boundary layer at 


{x, y) = 


Cr - Xo , Cr - Xo , 

7 7 H 

CT Lq^ CT 


(28) 


The tangent and normal unit vectors to the boundary layer are respectively given by 
t = {lx,ly)/\l\ and n = {ly, —lx)/\l\, where \l\ = + ly. Therefore, the normal and 

tangent components at (x, y) are calculated as follows; 



(> 0), 

(29) 

6 ^^ (^-^a + ly 


(30) 

(<0), 

(31) 

4“^ - d) 

(<0), 

(32) 


where := ipoa^. 

Next, we consider the down-stream side, on which the magnetic field lines are again 
approximated by straight lines as shown in Fig. Since the displacement map is area¬ 
preserving, the same field line that passes through (x, y) is found to be 


y = 


f CT 


U 


I 


-X a — Xq 


(33) 


by equating the areas of the two blue triangles in Fig. Using the fact that the value of xjj 
is again '0’^°^(xo, 2/o) on this held line, a straightforward calculation results in 


'ip^'^\x,y,t) [^y-T^ + 


.-d 


y 2/^ 


a 


T: 


21 


ly 2 


a 


H 


L 


2L 


-xy 


a X 


-ax 


L 


+ "00 - 


a 


(34) 


The held line (33) also moves in time because of the time dependence of lx,ly and a. By 
imposing the boundary condition 0(0, y,t) = 0 on the y axis, the associated incompressible 
how can be determined uniquely as 


(t)^^\x,y) = 


2(3 x^ (3 


a 


p 




where (3 = Ix/ly We thus obtain, at {x,y), 


bf I ? + 


‘H 




21 


: - q- \ ^ 
lx V i/r 


Xd) 


v: 


It 


- ||s J ^ 

(j ^ rn 


(35) 

(36) 

(37) 
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and confirm that bn^ = and are indeed satisfied. 

The speed Vn for movement of the bonndary layer at {x,y) is calculated by using the 
angle 9 between the boundary layer and the y axis (tan^ = Ix/ly), 


Vn = = i/?; 


(38) 


Now, we are ready to impose the matching conditions on these up- and down-stream 


solutions. It is interesting to note that the matching condition (22) is already satished 
because we have taken the continuities of "0 and 0 into account in the above construction. 


The matching condition (21) at (f, y) gives 


- d = -Tj^^XoC, (39) 

where C, '.= Sy — {3 = {a — lx)/ly This condition must be satisfied for all points (x, y) on the 
boundary layer (that is, for all xq G [lx, o']), which requires both (a = a( and ( = to 

be satisfied. The former gives a constant of motion. 


const., (40) 

(41) 

Although the constant lyo{> 0) is still unknown unless lx and ly are specihed, the displace¬ 
ment a turns out to grow explosively due to the presence of the X-shaped boundary layer. 
These parameters lx, ly and lyo will be determined later when this inner solution is matched 
with the external solution and the global energy balance is taken into account. 


C CT-lx^ 


and the latter gives an evolution equation. 


a = 


lyO^H 


a 


C. External solution 

Now consider (iii), the external solution of Fig. Even though we discuss the nonlinear 
phase, the displacement a (or the island half-width) must be small a Lx/4: as well as the 
growth rate a/a in comparison with the equilibrium space-time scale. Therefore, 

we expect the external solution to be similar to the well-known eigenfunction of the linear 
tearing mode. This treatment for the external solution is commonly used in Rutherford’s 
theory [23l [2l], while we introduce the arbitrary wavelength Ly{< Ly) of the linear tearing 
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mode in this work. Namely, the displacement in the x direction is given by ^{x,y) = 
e^{x) cos{27iy/Ay) for y G [0,Ay/4] where the eigenfunction, 



cos 


- (L^/Ay)'^ {a\x 
sin ax 



(42) 


is normalized so as to satisfy ^(La;/4) = —1 and ^{—Lx/4:) = 1. Taylor expansion of ^ around 
X = 0 on the positive side (x > 0) gives 


^^(x) 




{L^/Ay)‘^cos 




+ 0(x), 


(43) 


where the tearing index A' for the wavelength Ay is 


A' = 2aJl — {Lx/Ay)'^ taxi 


TT 


1 - {Lx/Ay)^- 


(44) 


Since the dependence of ^ on Ay is complicated, we again restrict the range of A' to the 
large A' regime. 


A' > a;, := 100 /L,, (or Ay/Lx > 2.5), 


(45) 


as we have done in the linear theory. Then, we can use the following approximations: 


^(x) ~ -1 

A 2 

A' ~ 16-|. 
Li 


A'x 


+ 0(x), 


(46) 

(47) 


The critical value A), = 100/La, is, of course, specihc to the equilibrium state 
This external solution is matched to the inner solution by 

2 


-a = e^(a) = -e 1 + 


A'a 


+ 0{ea), 


(48) 


which gives, by neglecting 0{ea), 


a =e- 




(49) 


Note that a is larger than e as illustrated in Fig. Since the displacement map is area¬ 
preserving, we determine ly by the relation. 


aly 



(50) 
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D. Energy balance 


In linear tearing mode theory, the released magnetic energy via reconnection is estimated 


by 




(51) 


in terms of the perturbed flux function ip = —{dip^^^/dx)^ = where 2a(<^ L^) is the 

width of the boundary layer at a; = 0 (see Appendix A of Ref. H). This is true if the 
eigenfunction ip{x) is smoothed out and flattened within the layer [—a, a] by some sort of 
nonideal MHD effects. 

For the nonlinear phase in question, we simply replace a by a (and Ly by A^^) because 
the flux function ip is flattened within [—cr, a] by the formation of a magnetic island. This 
idea is similar to the hnite-amplitude generalization of A' which is made by White et al. [21] 
for the purpose of introducing a saturation phase to Rutherford’s theory. In either case, the 


island width 2cr grows as far as ipdip/dx\%^_'P > 0. Using the Taylor expansion (46) and the 


relation (50), we obtain 


2 7 \x=o 


dw = -e^^ip^y = -Alyrd^a^. 
2^ dx ^ ^ 


(52) 


We remark that the magnetic energy in the area [0,cr] x [0,/y] at the equilibrium state 
{t = 0) is also of the order of lyTp^'^a^. Since this area is mapped to the internal region of 
the magnetic island after the displacement, we can expect a corresponding decrease in the 


total magnetic energy, which agrees with the estimation (52). 

In order to satisfy the energy conservation 6K + dW = 0, the kinetic energy is required to 
satisfy 6K cx lya^. To be concise, let us assume ly = const, a priori because this assumption 
turns out to yield the desired scaling 6K cx as follows. 

Since the kinetic energy is mostly concentrated on the down-stream side due to the outflow 
from the X-shaped vortex layer, we use (p^'^'> in (35) and the orderings ly 3> h and ly 7^ a to 
estimate the kinetic energy in the down-stream region as 


dy 


'ky 


_ 

“ 22 /’ 


( 53 ) 
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where we have neglected in the last expression since Ix/cr = 1 — ly/lyo is now constant 

and less than unity. The same estimate is more easily obtained as follows. Consider the flow 
passing through the box [0, l^] x [0, ly]. Since the inflow velocity into the box is at most &, 
the outflow velocity, say Vy, is roughly determined by the incompressibility condition, 


&ly '^y^X') 


(54) 


where S> d owing to ly S> lx- The kinetic energy density Vy/2 multiplied by the area 


lxly/2 reproduces the same estimate as (53). 

In fact, the outflow also exists over the area [0,4] x [ly,Ly/2] in Fig. and there are 
eight equivalent areas in the whole domain according to the parity. Therefore, a plausible 
estimate of the total kinetic energy is 


2 r 


5A- = 8 | 


I 7 f ^ 

2 ^ " V 2 




= V 


(h 


n-2FiJ 

''v) n ^2 1 

'-yO' H 


k 


ly/lyO 


(55) 


where the evolution equation (41) for a has been used. Since this 6K is proportional to 
we can impose the energy conservation law 6K + 5W = 0, which determines lyo with respect 

to ly ^ 


l + ,/ 2 ^ 


hfi _ _ 

ly 2 


(56) 


Given this lyo, the rate of decrease in the magnetic energy 

dt{6W) 3a 


5W 


lyO^H 


(57) 


indicates that the shorter the length ly, the faster the magnetic energy decreases. Thus, the 
local reconnection (i.e., the local X-shape) develops faster than the global one. 

However, there is a lower bound for ly since this argument is based on the assumption 
A', < A' < A', for which the approximation A' = 16Ay/L^ (and A' = IQLy/L^) is valid. 


Using (49) and (50) with this approximation, ly can be regarded as a function of e and A', 


L 


\/Lx A' 


=■■ F{e,A'), 


( 58 ) 


\X + 


and hence ly/Lx should he between F(e, A),) and F(e, A') as shown in Fig. 10 
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FIG. 10: Range of ly corresponding to A'^ < A' < A' for the case of ky/a = L^jLy = 0.171. 


E. Scaling of the explosive growth 


By rewriting (58) as 




(59) 


this relation is found to have two kinds of scaling depending on whether its right hand side 
is much smaller or larger than unity, 

First, when lyjL^ ^ the relation reduces to 


8 

= —— -C 1. 
eA' 32/2e 


(60) 


Since a ~ e in this case, we obtain the same explosive growth as (41), 

£2 
e = 


TnlyO ’ 


(61) 


in terms of the displacement e at x = L^/A. We refer this scaling as kink-type because A' 
is so large that the external solution is similar to the kink mode (cr ~ e, ~ ^y/^)- 
On the other hand, when ly/L^ ^ fhe relation (59) reduces to 























By noting that 


e 

(T ~ - 

2 V eA' 




1/4 




( 63 ) 


the explosive growth (41) becomes 


^2 


(64) 


We refer this scaling as tearing-type because A' is so small that the external solution is 
similar to the tearing mode (a > e, ly < Aj,/4). 

The boundary line ly/L^ = between the kink-type and tearing-type regimes is 


also drawn in Fig. 10 Since the magnetic energy is released more effectively for the smaller 
ly, the fastest reconnection occurs near the lower bound ly/ = 0.4 ~ 0.6. Figure 


11 


shows 


that the tearing-type scaling (64) for ly/ = 0.4, 0.5, and 0.6 agrees well with the simulation 


results. For comparison, we also draw the kink-type scaling (61) with ly = lAL^ — Ly/A 
as a global reconnection model, which is indeed slower than the simulation results. We can 
conhrm that the stream lines in Fig. [^a) are more like the tearing-type {a > e, ly < Ay/4). 
Although the current layers in Fig. |^c) are actually curved, they are locally regarded as 
straight lines around the origin as in Fig. [^and seem to have the length ly/ — 0.5. Note 
that the same ly/L^ — 0.5 is also observed in Fig. |^b) and (c) since this ly/ is determined 
independently of Ly and de = ps as shown in 


V. SUMMARY 

We have investigated the nonlinear evolution of a collisionless tearing mode that can grow 
explosively with the formation of an X-shaped current-vortex layer due to the coexistence 
of electron inertia and temperature effects, where we have assumed de = Ps for simplicity. 

For the equilibrium state given in ([^ and the wavenumber ky = 2tt/L y, the tearing mode 
is linearly unstable when the tearing index A' (which is a function of Ly/L^) is positive. The 
simulation results show that explosive growth occurs when A' > 0.65/ps. More specifically, 
the amplitude e of the displacement at x = LL^/A exceeds ps and then grows explosively; 
e oc €"■, n > 1. By observing this explosive phase in detail for de = Ps < 0.01, we hnd that 
the X-shaped layer widens locally around the reconnection point and its length scale (~ 2ly) 
seems to be unrelated to the wavelength Ly (and A' as well) of the linear eigenmode. 
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FIG. 11: Logarithmic plot of the displacement e versus its time-derivative e for de = Ps = O.OO 5 L 3 ;, 


ky/a = 0.171. Dotted lines are theoretically derived from the kink-type scaling (61) with lyjLx = 


1.4 and the tearing-type scaling (64) with ly/Lx = 0.4,0.5,0.6 


To explain this locally enhanced reconnection, we have developed a theoretical model in 
which the magnetic flux 'ijj is assumed to be conserved (like ideal MHD) except within the thin 
X-shaped layer. Namely, the two-fluid conservation laws (|^, (|^ are invoked only within the 
layer to obtain the matching conditions across it. The external solution is approximated by 
a linear tearing eigenmode that has a shorter wavenumber Ay than Ly and a smaller tearing 
index A' than A'. We have restricted our consideration to the range A(,(= lOO/La,) < A' < 
A' (or 2.5Lx < Aj, < Ly), in which a simple expression A' = 16A^/L^ holds and the length 
of the local X-shape (~ 2ly) is related to A' by (58). As shown in Fig. 10, we have found 


that there are two kinds of scaling depending on whether the external solution is kink-type 
or tearing-type. The faster reconnection is theoretically predicted at the shorter ly, namely, 
at the lower bound of this range, ly ~ 0.5Lx, A' ~ A(, and Ay ~ 2.bLx, which belongs to the 
tearing-type regime. The simulation results indeed agree with the tearing-type scaling with 
the explosive growth rate e oc and they corroborate other properties predicted by this 
local reconnection model. 

In comparison with the classical Petschek reconnection model [2] in which the X-shaped 
boundary layer is composed of stationary slow-mode shocks, our model suggests that the 
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X-shaped current-vortex layer is kinematically generated by ideal, incompressible and ac¬ 
celerated fluid motion in accordance with the two-fluid conservation laws and the energy 


conservation. The explosive growth rate (64) with ly ~ O.bL^ is moreover independent of 
the microscopic scale dg = Ps and hence reaches the Alfven speed e ~ Ll / {Tulyo) at the fully 
reconnected stage e ~ L^/A. This is faster than the explosive growth e ~ kyd]J‘^e ^/‘^/th that 
is caused by the Y-shaped layer in the presence of only electron inertia (see our previous 
work j25]h 

The two-field equations ([^ and ([^ can be derived from gyrokinetic and gyro-fluid equa¬ 
tions by taking the fluid moments and then neglecting the ion pressure and electron and 
ion gyroradii [HI HE] . This fact suggests our present results are a barebones model for fast 
reconnection, but further generalizations including the case dg 7 ^ Ps are suggested for future 
work. The existence of more than one microscopic scale gives rise to nested boundary layers, 
as already known from the linear analysis. The nonlinear evolution of such nested boundary 
layers would be more complicated than that of the single boundary layer {dg = ps) we have 
discussed. Nevertheless, if the outermost layer is sufficiently thinner than the island width 
and the energy balance is dominated by ideal MHD, we expect a similar X-shaped layer 
and explosive growth, since it is unlikely that any other structure can exist that is more 
efficient for releasing magnetic energy. Unfortunately, present computational resources are 
not enough to observe the explosive phase for a sufficiently long period in the presence of 
the nested bounded layers. For example, when dg ps, linear analysis indicates that the in¬ 
nermost layer width ~ is even narrower than dg and demands more computational 

grids. Further advancements in computational performance and technique will be essential 
for studies of explosive reconnections in more general collisionless plasma models. 
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